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ABSTRACT 

The Standing Accretion Shock Instability (SASI) is commonly believed to be responsible for large 
amplitude dipolar oscillations of the stalled shock during core collapse, potentially leading to an 
asymmetric supernovae explosion. The degree of asymmetry depends on the amplitude of SASI, but 
the nonlinear saturation mechanism has never been elucidated. We investigate the role of parasitic 
instabilities as a possible cause of nonlinear SASI saturation. As the shock oscillations create both 
vorticity and entropy gradients, we show that both Kelvin-Helmholtz and Ray leigh- Taylor types of 
instabilities are able to grow on a SASI mode if its amplitude is large enough. We obtain simple 
estimates of their growth rates, taking into account the effects of advection and entropy stratification. 
In the context of the advective-acoustic cycle, we use numerical simulations to demonstrate how the 
acoustic feedback can be decreased if a parasitic instability distorts the advected structure. The 
amplitud e of the shock deformation is es timated analytically in this scenario. When applied to the 
set up of lFernandez fc ThompsonI ()2009a| ). this saturation mechanism is able to explain the dramatic 
decrease of the SASI power when both the nuclear dissociation energy and the cooling rate are varied. 
Our results open new perspectives for anticipating the effect, on the SASI amplitude, of the physical 
ingredients involved in the modeling of the collapsing star. 

Subject headings: hydrodynamics — instabilities — shock waves — supernovae: general 



1. INTRODUCTION 

Des pite decades of a ctive research ([Colgate fc White! 
[19661: iBethe fc WilsonI [l985l ). the core collapse super- 
novae mechanism remains elusive. The failure of the 
most sophisticated ID models to explode the majority of 
massive progenitors (Licbendorfcr ct al. 2001) suggests 
that multidimensional effects are essential for a success- 
ful explosion. Understanding the hydrodynamical in- 
stabilities responsible for this symmetry breaking, and 
more specifically their nonlinear dynamics, is therefore 
required to understand the explosion mechanism. 

The region between the neutrinosphere and the 
shock deserves particular attention because two in- 
stabili ties take place the re: ne utrino-driven convec- 
tion (i Herant et all 119921 119941: [Bu rrows et~all 119951: 
iJanka fc Miilled 119961: iFoglizzo etal . 2006), and the 
newly disco vered Standing Accr e tion Shock Instabil- 
ity (SASI) (IBlondin et alT [2?M [Ohnishi et al.l ^2XM 
iFoghzzo et al.l 120071 : IScheck et all [200^ 2D simu- 
lations suggest that the complex fluid motions trig- 
gered by these instabilities could lead to a success- 
ful explosion either by helping the classical neutrino- 
driven mechanism (Buras ct al. 2006.; Marek & Janl^ 
[20091: iMurphv fc Burrowsl 120081 ) or by a new mecha- 
nism based on the emis sion of acoust ic wav es from 
the proto-n eutron star (iBurrows et al.l ([20061 . I2007D . 
see however [Weinberg fc QuataertI ([2008f )). The large 
scale {I = 1 — 2) induced asymmetry could also ex- 
plain the h igh kick velocities of ne wly formed neu- 
tron stars ([Scheck et all 12004 I2006D and may affect 
significantly their spin ([Blondin fc Mezzacappal 120071 : 
lYamasaki fc Foghzzo|[2008lV 

The linear phase of the two instabilities has 
been 



describe d in de ta ils b y 'Foglizz o et al. 
[Blondin fc Mezzacappal ([20061 ): .Foglizzo et al. 



(2006): 



20071) : 



lYamasaki fc Yamadal ([2007| ): iFernandez fc ThompsonI 
([2009a). Neutrino-driven convective modes with a large 
angular scale can be stabilized by a fast advection of 
matter through the gain region, whereas SASI is always 
dominated by large scale modes. This linear argument 
favors SASI as the cause of the prominent 1 — 1 — 2 shock 
oscillations observed in the simulations. However, a the- 
oretical understanding of the nonlinear development and 
saturation of SASI is still missing. This was highlighted 
by the unex pected dramatic decrease of the SASI power 
observed by IFernandez fc ThompsonI ([2009a|) when both 
the nuclear dissociation energy and the cooling rate are 
varied. To shed more light on this issue, this paper pro- 
poses a predictive saturation mechanism for SASI. This 
is a first step toward understanding how the amplitude 
of SASI depends on the physical ingredients of the model 
(nuclear dissociation, equation of state, heating rate, ro- 
tation, magnetic fields). 

We propose that the saturation takes place when a par- 
asitic instability (also called secondary instability) grows 
on the dominant SASI mode. These parasites feed upon 
its energy and destroy its coherence, leading to a turbu- 
lent flow and the saturation of the growing SASI mode. 
This saturation mechanism by parasitic instabilities is 
si milar to the one proposed for the saturation of t he MR I 
bv lGoodman fc Xul (|l99l : IPessah fc GoodmalU (|200^. 

Two types of instabilities are considered as potential 
parasites for the SASI mode: the Kelvin-Helmholtz insta- 
bility (hereafter KHi) grows near a maximum of vorticity, 
and the Rayleigh- Taylor instability (hereafter RTi) grows 
on negative entropy gradients. In either case the cause 
of the instability lies within the SASI mode, and there- 
fore the growth rate of these parasites increases with the 
amplitude of the shock oscillations. The parasites are 
thus able to affect significantly the dynamics of the flow 
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only if the amplitude of SASI oscillations exceeds a cer- 
tain threshold. The main objective of this work is to 
estimate this threshold. 

Understanding the saturation mechanism of SASI re- 
quires in principle some understanding of the mecha- 
nism underlying its linear growth. The saturation sce- 
nario we propose is not restricted to a particular mech- 
anism for the growth of SASI, but it can be under- 
stood more preci sely in th e fram e work of the advective- 
acoustic cycle (Foglizzo' '2002'; 'Blondin et al.' '20031; 
Ohni shi et al. 2006; Foglizzo et al. 2007; Sc hcck c t all 
20081 : (Fernandez fc Thompsonll2009at lFogliizoll2009h . In 
this mechanism, a shock deformation creates an entropy- 
vorticity wave, whose downward advection generates an 
acoustic feedback. This acoustic wave further deforms 
the shock, thus closing the unstable cycle. If the coher- 
ence of the advected wave were broken by a parasitic 
instability before it created the acoustic feedback, this 
cycle would be stabilized and SASI saturated. 

In Sect. 2, we explain our method and approxima- 
tions. Sections 3 and 4 are devoted to the linear study 
of the KHi and RTi in our setup. In Sect. 5 we study 
how the development of the parasites breaks the coher- 
ence of a SASI mode and decreases of the acoustic feed- 
back. In Sect. 6 we ap p ly our results to the setup of 
[Fernandez fc ThompsonI (|2009aD and compare our esti- 
mates with the results of their simulations. Finally, our 
results are discussed in Sect. 7 and summarized in Sect. 8. 

2. METHOD 

2.1. Estimating the growth rate of the secondary 
instabilities 

We estimate the local growth rate of the parasitic in- 
stability by using a simplified description of the linear 
SASI mode that keeps only the features that are essen- 
tial for the physics of the instability. 

A SASI mode has a complex structure of entropy, 
vorticity and pressure perturba tions which can be 
computed by a linear an alysis (|Foglizzo et al.l 120071 : 
lYamasaki fc Yamadal 120071) . We focus on the advected 
structure of the SASI mode, where the parasitic insta- 
bilities operate. As we seek a local description of the 
parasites, we assume a planar geometry. Denoting by 
z the vertical direction, and x the transverse direction, 
the structure of the advected entropy /vorticity wave (de- 
noted as S, w respectively) is approximated by a sinusoid 
with a vertical wave number Kz = i^/vq, where w is the 
SASI frequency and vq is the flow velocity: 

S{z,t) = So{z) + ASxcos{ujt-Kzz) (1) 
w{z,t) ^ Aw X cos{u}t ~ Kzz), (2) 

where Sq is the entropy profile of the radial flow, and 
AS", Aw are the amplitudes of entropy and vorticity 
perturbations associated with the SASI mode. The ad- 
vected entropy/ vorticity wave in a SASI mode is actu- 
ally tilted with respect to the horizontal direction, but 
we neglect this tilt for the sake of simplicity {K^ is 
a factor > 5 larger than the horizontal wave number 
K^jy + 1) /r in a typical SASI mode (jFoglizzo et al.l 
l2007f )). In what follows we omit the index z and simply 
note K = uj/vq the wave number of advected SASI per- 
turbations. The horizontal wave number of each para- 
sitic instability is denoted as k. The notation Aq refers 



to a quantity (A) in the stationary flow, unperturbed 
by SASI. Its perturbation by the SASI mode is denoted 
as AA, and 6A refers to a parasitic perturbation. The 
growth rates of the KHi and RTi are directly related to 
the amplitude of the SASI mode through its profile of 
vorticity and entropy gradient. Non adiabatic cooling is 
taken into account in the shape of SASI eigenf unctions, 
but is neglected in the dynamical evolution of the para- 
sites. We choose the usual adiabatic index of a relativistic 
gas 7 — 4/3. 

We estimate the stability of the fiow described by 
Eq. ([1]) and ^ in three steps, first by neglecting the 
background entropy gradient and global advection vq, 
second by assessing the effects of a uniform entropy gra- 
dient VS'o, and third by taking into account advection. 

In the first two steps, a standard linear mode analysis 
can be used because the flow is stationary. The equations 
determinin g the evolution of the pert urbations (e.g. Ap- 
pendix A of lFoglizzo fc Ruffertl (fToM ll are solved numer- 
ically. A simple fitting formula for the maximum growth 
rate is proposed, as a function of the SASI amplitude and 
the background stratification. 

The third step is related to the concept of gl obal ver- 
sus local instability (jHuerre fc Monkewitzill990l ) , but the 
gradients in the direction of advection preclude the use of 
standard analytical techniques. If the fluid is advected 
too fast, the instability may be able to grow in a la- 
grangian way but would actually decay at a fixed radius 
as the perturbations are advected away. We use numeri- 
cal simulations to measure the propagation speed of the 
parasitic instability in the z-direction, by perturbing the 
SASI mode over a limited region. Adjusting the mea- 
sured speeds with a physical but approximate description 
of the propagation, we obtain a simple analytical esti- 
mate of the growth rate cTparasito of each parasitic insta- 
bility, taking into account the SASI amplitude, the back- 
ground stratification, and the advection speed (Eqs. 1121 
and nil). 

2.2. Estimating the saturation amplitude 

The growth rate CTparasitc is an increasing function of 
the SASI amplitude. A parasitic instability can affect the 
dynamics of SASI if its amplitude 6A becomes compa- 
rable to the SASI amplitude AA. The growth of the 
ratio SA/AA requires cTparasito > Csasi- We use this 
criterion to estimate the s aturation amplitude of SASI. 
iPessah fc GoodmanI ([20091 ) use a similar criterion to es- 
timate the saturation amplitude of the MRI. 

The criterion CTparasitc = CTsasi defines the minimum 
amplitude Aylmin(r) of SASI above which parasites can 
compete with SASI at a given radius r, despite advec- 
tion and cooling. The parasitic instabilities can alter 
the growth of SASI only if their growth takes place in 
a region which is vital to the mechanism of SASI. For 
example, if the mechanism of SASI is interpreted as an 
advective-acoustic cycle, this cycle is most sensitive to 
the region between the shock and the deceleration region 
where most of the acoustic feedback is produced. Fortu- 
nately, as will be shown in Sect. 6, the local saturation 
amplitude AA^in{r) displays a broad minimum around 
the radius (r* -|- rsh)/2, which defines a global saturation 
amplitude AAmin without much sensitivity on the details 
of the SASI mechanism. 
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2.3. Limitations 

By focussing on the growth of parasitic instabilities, 
we ignore other nonhnear processes which could play 
a role in saturating the amplitude of SASI. Among 
them, the steepening of acousti c waves into shocks (e.g. 
[Fernandez fc ThompsonI ()2009b[ )). the decoherence of the 
mode due to the finite displacement or velocity of the 
shock, or the exchange of energy by resonant mode cou- 
pling, could be important in some parameter range. 

When compared to published simulations (Sect. 6), our 
estimate of the saturation amplitude of SASI based on 
parasitic instabilities is encouraging in view of the many 
simplifications inherent to our method: our description 
of the parasitic growth neglects the spherical geometry 
of the flow, neglects the tilt of the SASI wave with re- 
spect to the horizontal direction, and assumes an adia- 
batic evolution of the parasites. We further assume that 
SASI is dominated by a single SASI mode of finite am- 
plitude, which we describe using a linear approximation. 
Steepened advected waves, induced by steepened acous- 
tic waves reaching the shock, could affect the growth of 
parasitic instabilities by introducing larger vorticity and 
sharper entropy gradients. We also neglect the produc- 
tion of entropy by acoustic waves steepening into shocks, 
and the creation of vorticity by the baroclinic interaction 
of entropy gradients with pressure waves. 

3. THE KELVIN-HELMHOLTZ INSTABILITY (KHI) 

3.1. The KHi in a sinusoidal velocity profile 

The KHi feeds on the kinetic energy available in shear 
flows. A necessary condition for its growth is the pres- 
ence of a maximum in the absolute value of vorticity 
(jPrazin fc Reidlll98H ). The linear growth of SASI cre- 
ates a sinusoidal velocity profile 

v.,{z)= Av sin (Kz), (3) 

with two suc h maxima per wavelength. 

iHevvaerts fc PriestI (|1983f) have demonstrated that 
this profile is indeed unstable to perturbations with a 
small horizontal wave number k < K, with a growth 
rate a < kAv. 

In order to estimate the maximum growth rate and the 
corresponding wavelength, we complement their analyti- 
cal study by a numerical mode calculation where the flow 
profile is described by Eq. ([3]), with periodic boundary 
conditions in z = and z — 1. The results (Fig. [1]) show 
good qualitative ag reement with the expectation from 
IHevvaerts fc PriestI (|1983f ): the maximum growth rate 
(7 = 0.27i4rAu is reached for a wave number k = 0.58^, 
and the KHi is stable for k > K. The maximum KHi 
growth rate is thus a fraction of the max i mum vortic- 
ity in the flow (see also iFoglizzo fc RuffertI ()1999l )). The 
Mach number Ai has little effect on the instability unless 
it approaches unity (we used Ai = 0.2 in our mode calcu- 
lation, with results almost identical to the incompressible 
case). 

3.2. Effect of stratification on the KHi 

The buoyancy force in a stably stratified atmosphere is 
able to stabilize the KHi if Ri > 1/4, where the Richard- 
son number Ri characterizes the relative strengths of 
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Fig. 1. — Growth rate a of the KHi as a function of the transverse 
wave number k. The dashed hne is the eigenvalue for a sinusoidal 
velocity profile and periodic boundary conditions. The full lines 
show the effect of stratification with the Richardson number Ri 
ranging from to 0.2. The triangles and circles are measured 
from simulations where perturbations are localized (triangles) or 
extended (circles), a is normalized by the maximum vorticity ui, 
and fc is normalized by the vertical wavenumber K of the SASI 
mode. 

buoyancy and shear (e.g. iChandrasekhaii (|1961t )): 



where w = V x w is the vorticity. In the absence of a 
composition gradient, the Brunt- Vaisala frequency iV is 
defined by: 

iV^ = _l^gV5. (5) 
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The entropy S is here measured in dimensionless units 
5^1og(P/p^)/(7-l). 

In order to characterize this stabilization quantita- 
tively in our specific geometry, we compute the eigen- 
modes of the KHi in a sinusoidal horizontal velocity field, 
embedded in a stable entropy gradient and a gravity field. 
The sinusoidal velocity profile, the background entropy 
gradient and gravity profiles have a constant amplitude 
in a limited region of space and are smoothly connected 
to zero outside this region. This region has a size of 3 
vertical wavelengths. If far enough from the gradients, 
the choice of boundary condition (reflective or leaking 
condition) does not change the resulting mode. Results 
are shown in Fig. [TJ The slight difference between the 
curve where the entropy gradient is zero (Ri = 0) and 
the periodic KHi is due to the limited size of the entropy 
gradient and gravity profiles. 

The maximum growth rate decreases linearly to zero 
as Ri is increased, while the wave number at maximum 
varies from k ~ 0.6i^ to fc ^ 0.8iir (Fig.[T]). The marginal 
stability is reached for a critical value of the Richardson 
number Ri = 0.24 in close agreement with the expected 
value Ri = 1 /4. 

3.3. Effect of advection on the KHi 

In the context of core collapse, the SASI vorticity wave 
is advected toward the neutron star. Using a physical 
argument, we first give a naive estimate of the speed at 
which an unstable perturbation can propagate against 
the stream, and evaluate the reduced growth rate in the 
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presence of advection. Numerical simulations are used in 
order to obtain more accurate estimates. 

3.3.1. Phystcal argument 

The simple st illu stration of the KHi (e.g. 
iDrazin fc ReidI ()1981l )) considers an incompressible 
fluid with a discontinuity of horizontal velocity: 
Vx = —Av for z < 0, and Vx = Aw for z > 0. Unstable 
modes exist for any horizontal wave number with a 
growth rate (Tq = kxAv, and their structure on both 
sides of the discontinuity is a decreasing exponential: 
SA oc e~*''==l^l. Viewed in a frame moving with a vertical 
velocity Vz, the time dependence of the perturbation 
is: SA oc ei<^o-k^v^)t given height z above the 

discontinuity. A condition for the perturbation to be 
growing at a given radius can be deduced: ctq > kxVz, 
which may be interpreted as a propagation speed fprop 
of the KHi equal to: 

Wprop = TT- (6) 

The local growth rate in the frame moving with respect 
to the discontinuity is thus decreased as follows: 

a^ao(l^^). (7) 

Of course this physical argument is too simple to be 
directly applicable to the case of a sinusoidal shear wave 
advected downward, where each unstable shear layer is 
followed by another one. Could the propagation be faster 
if these adjacent shear layers cooperate ? 

3.3.2. Numerical simulations of the KHi 

Using the code RAMSES (|Tevssieri [2Ciol 
iFromang et al.l I2OO60 we performed numerical sim- 
ulations of a sinusoidal velocity profile described by 
Eq. © (with = 27r, Aw = 1, c = 5). The compu- 
tational domain was a box —8 <z<8, 0<a;<8, 
with periodic boundary conditions. Perturbations of the 
vertical velocity, localized between z — —0.5 and z — 0.5 
(one wavelength), have been added to the stationary 
flow, with an amplitude of lO^^Aw. The runs presented 
here have a resolution of 1024 x 512, but we checked 
that the results do not depend on the resolution. 

The different phases of the KHi growth are illustrated 
by three snapshots in Fig. [2] The wavelength ~ 1.8 of the 
dominant feature agrees with the expected wavelength of 
the fastest mode in Fig. [T](fc/iir ~ 0.55). As a consistency 
check, the growth rates have been measured for different 
wave numbers and compared successfully with the linear 
analysis in Fig. [J]. 

In the two bottom plots in Fig.[2]the KHi mode has al- 
ready propagated downward and upward from the initial 

^ The growth rates (triangles) are shghtly smaller than the linear 
values with periodic boundary condition, hut are close to those of 
the modes whose spatial extent is restricted to three wavelengths 
(Sect. 3.2). Note that they are not expected to match exactly, but 
the fact that they have similar values suggests a common physical 
origin of the reduced growth rate, namely the restricted extent of 
the region where the KHi grows. For comparison we also ran a 
simulation where the whole flow was perturbed instead of the re- 
gion z = [—0.5,0.5]. The growth rates (circles) are in very good 
agreement with the modes computed using periodic boundary con- 
ditions (dashed line). 
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Fig. 2. — Different stages in the evolution of the KHi on a sinu- 
soidal velocity profile. The upper plot is the initial condition: a 
sinusoidal transverse velocity wave with random perturbations lo- 
calized between —0.5 < z < 0.5 (one wavelength of the stationary 
flow). The middle plot shows the time when the KHi just reached 
a nonlinear amplitude, the mode structure is still clear. Finally 
the bottom plot shows a more developed nonlinear stage of the 
instability. In all plots, the grayscale represents the vorticity. 
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Fig. 3. — Propagation of the KHi along the z direction. The 
three curves represent different modes with a number of horizontal 
wavelengths ranging from 2 (upper curve) to 6 (lower curve). The 
plotted quantity, 2+ (t), corresponds to the spatial extent in which 
the KHi has reached an amplitude of Sv = W'^Av. The curves 
were shifted in time so that they start to deviate from at t = 0. 
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Fig. 4. — Upper plot: Propagation speed Vprop of the KHi as a 
function of tlie wave number k, measured in the numerical simu- 
lations (triangles). The full line shows the quantity 2.8<j/k, where 
cr is the growth rate computed in the linear analysis of Sect. 3.1. 
The factor 2.8 has been adjusted in order to match the results of 
the simulations. Bottom plot: maximum growth rate (full line, 
left axis) and associated wave number (dashed line, right axis) as 
a function of the advection velocity. The thin straight line with 
arrows is a linear approximation of the growth rate (Eq. ^) . 

perturbations. For a quantitative study of the propaga- 
tion speed, we measured the spatial range [— in 
which the KHi has reached an arbitrary amphtude (say 
Sv = 10^"* Au), as a function of time for different wave 
numbers. As expected from 3.3.1, the large wavelength 
perturbations propagate faster than the short wavelength 
ones (Fig. El). 

The propagation velocity measured in the simulations 
is typically a factor 2.8 higher than the quantity cr/fc 
suggested in 3.3.1, but the dependence on k is similar 
(upper plot in Fig.|3|). We interpret this high propagation 
speed as the sign that adjacent shear layers do cooperate. 

Approximating the propagation speed by Wprop = 2.8 x 
a/k, the wavelength dependence of the growth rate for 
a given advection velocity is estimated using Eq. ([7]). 
The maximum growth rate is smaller, and obtained for 
a longer wavelength than without advection, because 
shorter wavelengths propagate more slowly (Fig.]?]). The 
maximum growth rate dKHmax decreases almost linearly 
with the advection velocity Vz- We approximate this 
curve linearly as follows: 



CKHe 



(8) 



where ctokh = 0.25Aw is the maximum KHi growth 
rate without advection, and Veso is an effective propa- 



gation speed which we estimate to be (in the absence of 
stratification) VcSo ^ 1-5 x At; ~ 1.5Aw/K. Note that 
I'prop in Eq. ^ depends on the wavelength k of the per- 
turbation, whereas VeS in Eq. ([8]) is independent of k, 
because (TKHmax is the growth rate maximized over all 
wavelengths. Vee can be interpreted as an average prop- 
agation speed of the modes growing fastest at different 
advection speeds. 

3.4. Analytical estimate of the KHi growth rate with 
both advection and stratification 

In the presence of stratification, Eq. ([8]) becomes : 



CTKH — CTKHstrat 1 



(9) 



where UKHstrat = o-QKH (1 — Ri/Rio) is the maximum 
growth rate of the KHi in the presence of stratification 
but in the absence of advection, and v^s is the effec- 
tive propagation velocity in the presence of stratification. 
Eq. ([6]) suggests that the propagation speed is propor- 
tional to the growth rate, giving : 

fKHstrat /,„\ 

VcS = ■^effO, (10) 

CTOKH 

Injecting Eq. (|10l) into Eq. ^ then gives the following 
formula : 



, , Ri Vz 



(11) 



where ctokh = 0.25Aw is the growth rate in the absence 
of stabilizing effect, Rio = 0.24, and Vcso — l-5Aw/K. 
Equation (fTTj) can be rewritten as a function of the SASI 
amplitude: 



CTKH = 0.25Aw - 1.04 



Nl 
Au 
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(12) 



In a SASI eigenmode, the vorticity Aw at a given po- 
sition in the fiow is directly proportional to the relative 
shock displacement Ar/rgh- The factor Wsasi{r) defined 
by Aw = WsasiAr/rsh depends on the radius, and is a 
result of the linear mode analysis of SASI. 

A local saturation amplitude of SASI due to the par- 
asitic growth of the KHi, at a given radius, is deduced 
from Eq. ([T2|) and the criterion ctkh = cTsasi : 



'^Kvz 



2(7s 



1/2 



rsh 



(13) 

The discussion of a global saturation amplitude is post- 
poned to Sect. 6. 

4. THE RAYLEIGH-TAYLOR INSTABILITY (RTI) 

4.1. Simple RTi in a sinusoidal entropy profile 

The RTi feeds on the potential energy available when 
a low entropy fiuid is sitting on top of a higher entropy 
one, in a gravitational acceleration g. Buoyancy is char- 
acterized by the Brunt- Vaisala frequency N (Eq. (O). 
The flow is unstable if N"^ is negative, and the typical 
growth rate ctrt of short wavelengths perturbations is: 



-N 



2\ 2 
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gVS 



(14) 
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Fig. 5. — Growth rate of the RTi as a function of the transverse 
wave number. The full lines are the result of the linear analysis of a 
sinusoidal entropy profile embedded in a stable background entropy 
gradient and a gravity field. The different curves show different 
values of the background entropy gradient ranging from VSo = 
(upper curve) to V5o = 0.8V A5 (lower curve). The triangles show 
the growth rate measured in the numerical simulations without 
background entropy gradient. 

In the sinusoidal entropy profile created by SASI, the 
sign of the entropy gradient VS* changes every half wave- 
length: the entropy wave is made of adjacent layers of 
stably stratified and Rayleigh- Taylor unstable fluid. The 
RTi is expected to grow fastest where the entropy gradi- 
ent is most negative. 

By computing the eigenmodes of a sinusoidal entropy 
profile with a vertical wavenumber K m a constant grav- 
ity field, we verified that the growth rate a ^ ctrt at 
short wavelength (Fig.[5|), and a ~ O.TSctrt if fc ~ i^. In 
this calculation, as in Sect. 13. 2[ the entropy gradient and 
gravity profiles have a constant amplitude over a lim- 
ited region of space (three vertical wavelengths) , and are 
smoothly connected to zero outside this region. 

4.2. Effect of stratification on the RTi 

Just as KHi, the RTi can be stabilized by the presence 
of a stable entropy gradient in the stationary flow V5o. 
The flow is stable where the stationary gradient V^o is 
stronger than the SASI gradient VA5. The maximum 
growth rate is expected to scale as the Brunt- Vaisala fre- 
quency associated with the most negative entropy gradi- 
ent: 

^-gV{AS + So) \ 



CRT 
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(15) 



This was checked by adding a background positive en- 
tropy gradient to the previous mode analysis. The re- 
sulting growth rate, shown in Fig. [SI follows Eq. (|15l) at 
short wavelengths. At longer wavelength, the RTi is sta- 
bilized slightly faster. The cause may be that these larger 
scale modes, in addition to grow on less intense negative 
entropy gradient, are also sensitive to the higher positive 
entropy gradient (which are stable). 

4.3. Effect of advection on the RTi 
4.3.1. Phystcal argument 

The simple physical argument applied to the KHi in 
Sect. 3.3.1 can be adapted to the case of the RTi by con- 
sidering an incompressible fluid with a discontinuity of 
density with p = pQ — Ap for 2; < and p = po + Ap for 
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Fig. 6. — Different stages in the evolution of the RTi on a si- 
nusoidal entropy profile. The upper plot is the initial condition: 
a sinusoidal entropy wave embedded in a constant gravity field, 
with random perturbations of the vertical velocity localized be- 
tween —0.2 < z < 0.2. The middle plot shows the time when the 
RTi just reached a nonlinear amplitude, the mode structure is still 
clear. Finally the bottom plot shows a more developed nonlinear 
stage of the instability. In all plots, the grayscale represents the 
entropy. 

z > 0. Perturbations with a horizontal wavenumber kx 
are unstable with a gr owth rate ctq = {k^gAp/ p)^/"^ (e.g. 
iChandrasekhaH (|1961f )). The vertical structure of the 
RTi mode on both sides of the discontinuity is a decreas- 
ing exponential: 5 A oc e~'^="'^' due to the incompressible 
nature of the flow (k1 + =0). Using the same argu- 
ment as in Sect. 13. 3| the estimated propagation speed is 
Wprop = '^o/kx- In an entropy wave, each unstably strat- 
ified layer is followed by a stably stratified one. Unlike 
the instability of a vorticity wave, adjacent layers are not 
expected to cooperate. The vertical propagation of the 
RTi is expected to be accelerated by unstable layers and 
decelerated by stable ones. 

4.3.2. Numerical simulations of the RTi 

In order to measure the propagation speed of the RTi, 
we performed numerical simulations of a sinusoidal en- 
tropy profile S{z) described by: 



S{z) = ASsm{Kz) 



(16) 



In units of the vertical wavelength 2'k/K = 1, we chose a 
computational domain —8 < z < 8 and < a; < 8, with 
periodic boundary conditions in the x direction, and a 
resolution of 1024 x 512. The entropy oscillations de- 
scribed by Eq. (fTH)) and the gravity profile have been 
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Fig. 7. — Propagation of the RTi along the z direction. The 
three curves represent modes with a wave number ranging from 
k = Q.5K (upper curve) to A: = 2K (lower curve). The plotted 
quantity, (t), corresponds to the spatial extent in which the 
RTi has reached an amplitude of Sv = 10~*c. The curves have 
been shifted in time so that they start to deviate from at t = 0. 

restricted to — 5 < z < 5, the medium being uniform for 
1^1 > 5. At |z| = 8 we imposed a zero gradient boundary 
condition. The different phases of the RTi growth are 
illustrated by three snapshots in Fig. HI The RTi grows 
at all scales resolved by the grid. Numerical convergence 
has been checked for the modes studied in the following 
{k/K > 0.5), and the growth rates agree well with the 
linear analysis (Fig. [5]). 

Perturbations of the vertical velocity localized in a nar- 
row layer —0.2 < z < 0.2 have been added to the station- 
ary flow with an amplitude of 0.01% of the sound speed 
Co in order to follow their propagation. 

As expected in the previous Section, the propagation 
is less regular than for the KHi: it is slow in the sta- 
bly stratified regions and faster in the unstable ones 
(Fig. [71). It is however possible to measure a global 
propagation speed, which is remarkably close to the es- 
timate a/k deduced from the mode analysis (upper plot 
in Fig. |H]). The propagation speed is well approximated 
by ■yprop = 1.25cr/fc. 

Using Eq. ([7]), the RTi growth rate for a given advec- 
tion speed is expressed as a function of the wave num- 
ber. It is maximum at a longer wavelength than in the 
absence of advection, because shorter wavelengths prop- 
agate more slowly. We approximate linearly the decrease 
of the maximum RTi growth rate CTRXmax when the ad- 
vection velocity Vz increases (bottom plot in Fig. [8]): 



fRTn 



— CTORT 1 — 



WcffO 



(17) 



where ctort = 0.75AA^ is the RTi growth rate with- 
out advection and Veso is an effective propagation speed 
which we estimate to be ?;effo — 1-25AN/K. 

4.4. Analytical estimate of the RTi growth rate with 
both advection and stratification 

By a similar reasoning as in Sec. 3.4, the RTi growth 
in the presence of both advection and stratification is 




V /(N/K) 



Fig. 8. — Upper plot: propagation speed of the RTi as a function 
of its wave number k. The triangles are measured in the simu- 
lations. The black line is the quantity a/k estimated from the 
linear mode calculation. Bottom plot: maximum growth rate of 
the RTi as a function of the advection velocity (full thick line) and 
corresponding wave number (dashed line). The thin line with ar- 
rows ill ust rates the linear approximation of the growth rate used 
in Eq. lITTt . 



approximated using the results of Sects. |4?2] and 1473 
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-gV {AS + 5o) 



- 0.6Kvz. (18) 



The entropy gradient V (AS*) in a SASI mode is propor- 
tional to the relative shock displacement Ar/fgh- Defin- 
ing VS'sasi by V(AS') = VS'sasiAr/rsh, we use the cri- 
terion (jRT = CTsasi to obtain an explicit estimate of the 
local saturation amplitude of SASI due to the parasitic 
growth of the RTi: 



7 



Ar _ V^o ^ 
rsh VS'sasi 0.56 (7 - 1) gVSs, 



-I- OMsi^siVz) ■ 

(19) 



5. ACOUSTIC FEEDBACK IN THE PRESENCE OF 
PARASITIC INSTABILITIES 

The analytical estimates of the saturation ampli- 
tude obtained in Eq. ([T^ and ([TO]) from the criterion 
Cparasitc = Csasi are dircctly compared to the numeri- 
cal simulations of SASI in Sect. 6. Before that, we use 
the simplified toy-model of Sato et al. (2009) to evalu- 
ate the nonlinear effect of the parasitic instabilities on 
the ad vective- acoustic cycle. The distortion of the SASI 
mode by growing parasites, illustrated by the bottom 
plots in Fig. [2] and m is expected to induce a decrease in 
the acoustic feedback and stabilize the advective-acoustic 
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Fig. 9. — Efficiency of the acoustic feedback as a function of the 
amphtude of the advected wave. The dashed hne corresponds to 
an entropy vorticity wave (case i) , the full line to an entropy wave 
(case ii) and the dotted line to a vorticity wave (case iii) . The two 
vertical lines represent the estimate of the amplitude cutoff where 
the corresponding efficiency has decreased by 50%: dashed for the 
entropy-vorticity wave, full for the entropy wave. 

cycle responsible for the growth of SASI. 

The "problem 1" studied by iSato et alj (|2009D deals 
with the deceleration of an advected wave through an 
external potential, in a planar toy-model. This deceler- 
ation region of size A^v generates an acoustic feedback, 
measured at a distance Zmcas above it. The advected 
wave of amplitude eg is perturbed by a random noise 
acting as a seed for the parasitic instabilities. We choose 
Azv = 0.4, Ml = 5, cfjcl^^ = 0.75, ujt,,,/2tt = 2. 
The numerical technique based on a AUSMDV scheme 
is described in Sato et al. (2009). 

The effects of the KHi and RTi on the acoustic feed- 
back are studied together and separately by performing 
three sets of simulations: i) with the same mixture of 
entropy and vorti city as p roduced by a perturbed shock 
(Eq. (9-13) of Sat o et al.l ([2009)), ii) with the same en- 
tropy structure but no vorticity, iii) with the vorticity 
structure of (i) but no entropy. 

The pressure measured at Zmoas = 3 is Fourier trans- 
formed in the x-direction and in time, in order to esti- 
mate the part of the acoustic feedback which is coherent 
with the initial advected wave. This coherent feedback is 
responsible for the closure the advective-acoustic cycle. 

Above a certain amplitude threshold, the acoustic feed- 
back efficiency decreases from the value predicted by the 
linear analysis to a small fraction of this value (Fig. 
This threshold is measured as the amplitude 65 at which 
the acoustic feedback efficiency is 50% of its linear value. 
We find a value of eg = 0.68 for case (i) and £5 = 1.3 for 
case (ii). The decrease of the feedback is due to the de- 
velopment of parasitic instabilities propagating against 
the flow (Figs. [TU] and [T^ . The growth of the para- 
sites causes the advected wave to lose its coherence: the 
power remaining in its n = 1 component is decreased sig- 
nificantly in the region where the parasites have grown 
(Fig. [TT|) . The vertical structure of the feedback is also 
distorted. 

In Fig. ini the small amplitude of the acoustic feed- 
back in the case (in) decreases by 30% for es = 0.8, and 
increases again for eg = 1: this increase is due to the 
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Fig. 10. — Acoustic feedback from an entropy-vorticity wave (case 
i) in the presence of the KHi parasites. The left column represents 
a wave of linear amplitude (e^ = 0.1, below the cutoff), while the 
right column represents a wave of non-linear amplitude {es = 1, 
above the cutoff). The three rows show entropy (upper), vorticity 
(middle), and pressure (bottom) perturbations. The horizontal 
dashed lines represent the extent of the potential jump. The KHi 
is able to grow only on the non-linear wave (right column). 

pressure associated with the KHi, propagating against 
the fiow. 

Can the linear description of the KHi and RTi (Sect. 3 
and 4) predict the value of the threshold? We use the 
analytical estimates of the KHi and RTi growth rates 
(Eqs. (fT2|) and (fT8|) ) at marginal stability to estimate 
the threshold amplitude es{z) above which the KHi or 
RTi can grow despite the stabilizing effect of advection. 
The eg-threshold for neutral stability is expected to be 
a lower bound for the threshold measured in the simula- 
tions (Fig.Hni). 

The RTi can grow only in the region —0.2 < z < 0.2 
where gravity is significant, represented by the two hori- 
zontal dashed lines in Figs.[TU]to[T21 The local amplitude 
threshold in Fig. [13] displays a clear minimum around 
z ~ at a value of es — 0.78, which appears to be a 
good estimate of the amplitude above which the RTi can 
grow and damp the acoustic feedback. 

In the case (ii) where the upstream advected wave con- 
tains no vorticity, the growth of the KHi is subdominant 
and does not affect the acoustic feedback because it takes 
place below the region of deceleration (Fig. [13]). The RTi 



The saturation of SASI 



9 




1 2340 1 2J4 
X 



Fig. 11. — Same as Fig. [TO] but filtered to keep only the n = 1 
component of the horizontal structure. In the entropy and vorticity 
profile a clear decrease in amplitude is visible where the KHi has 
grown. In the right column, the pressure perturbations are slightly 
smaller and the coherence of the vertical structure is lost. 

is thus the dominant instabihty. The hnear threshold 
(eg = 0.78) deduced from Fig. [13] is about 40% smaller 
than the value of eg = 1.3 measured in the simulations. 

In the case (i) , the local amplitude threshold of the KHi 
measured in Fig. [13] is smaller than that of the RTi. The 
KHi should thus be the dominant parasitic instability 
with a smaller amplitude threshold than in case (ii) (eg e 
[0.36,0.65]). Indeed an instability develops upstream of 
the potential jump in case (i) with vortices which are 
typical of the KHi (Fig. [TUl) . For comparison, vortices are 
less prominent in case (ii) and no disturbance is growing 
upstream of the shock, consistent with the RTi (Fig.[T2]). 
As predicted the cutoff in case (i) is smaller than that of 
case (ii) (0.68 versus 1.3). The linear threshold (^ 0.5) 
is roughly 25% smaller than the simulated value of £5 = 
0.68. 

As a summary, wc find that the linear thresholds for 
marginal stability are consistently 25 — 40% smaller than 
the threshold at half efhciency measured in the simula- 
tions. Let us emphasize that in case (i), although the 
linear acoustic feedback is essentially generated by the 
entropy wave, vorticity plays an essential role in deter- 
mining the saturation threshold through the KHi. This 
illustrates the non-trivial interplay of vorticity and en- 
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Fig. 12. — Acoustic feedback from an entropy wave (case ii) in 
the presence of RTi parasites. The left column represents a wave of 
linear amplitude (eg = 0.1, below the cutoff), while the right col- 
umn represents a wave of non-linear amplitude (eg = 2, above the 
cutoff). The two rows show entropy (upper) and pressure (bottom) 
perturbations. The horizontal dashed lines represent the extent of 
the potential jump. The RTi is able to grow only on the non-linear 
wave (right column). 
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Fig. 13. — Amplitude eg corresponding to the marginal stability 
of seconda. r y inst abilities, as a function of z in the problem 1 of 
ISato et al.l RooSI ). Full lines show the KH i marginal stability: the 
entropy- vorticity wave (thick line, case i), the entropy wave (thin 
line, case ii). The dashed line shows the RTi marginal stability in 
both cases. The two vertical dotted lines show the vertical extent 
of the potential jump. 

tropy in the advective-acoustic cycle. 

6. COMPARISON WITH (MORE) REALISTIC 
SIMULATIONS: THE EFFECT OF NUCLEAR 
DISSOCIATION 

In order to test the parasitic scenario against results 
from numerical simulations, w e apply th e above esti- 
mates to the set up of iFernandez &: Thompson. (,2009al . 



10 



Ar/(r -r.) KH 




Fig. 14. — "Local saturation amplitude" as a function of the ra- 
dius for the KHi (upper plot) and RTi (bottom plot) instabilities. 
The thick lines show the fundamental mode for different dissocia- 
tion energies (upper curves correspond to lower energies) . The thin 
lines show higher harmonics with the dissociation energy e = 0.2ii^. 

where the energy loss at the shock (e in their notations) 
due to the dissociation of iron i s varied in a parameter- 
ized way from e = (the setup of lBlondin fc Mezzacappal 
(|2006[ )) to e = 0.25i;|. Here vg is the free fall velocity 
at the shock. The shock radius is kept constant by ad- 
justing the cooling function, which is varied by a factor 
up to 127 (the shock radius and the adiabatic index are 
''sh/'^star = 2.5 and 7 = 4/3). These simulations are a 
very good test for any saturation mechanism, because the 
saturation amplitude of SASI was found to be sensitive 
to the parameter e. 

By solving the radial structure of eigenmodes in the 
setup of Fernand ez &: ThompsonI p009a ). we calculate 
the parameters Wsasi and V5sasi and use Eqs. (IT^ and 
Uni) to estimate the "local saturation amplitude" of SASI 
oscillations above which the parasites grow faster than 
SASI at a given radius. As indicated by Fig. [TH the sat- 
uration amplitude of SASI associated with each parasitic 
instability decreases strongly when e increases. 

The global saturation amplitude can be estimated as 
the minimum of the local saturation amplitude, at least 
if this minimum is sufhciently broad and above the cou- 
pling radius. The curves in Fig.[T4]show a minimum at an 
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Fig. 15. — Roles of the advection, the entropy stratification, and 
the SASI growth rate in the "local saturation amplitude". This 
is illustrated with the fundamental SASI mode at e = and the 
RTi, but it is qualitatively the same if one considers other SASI 
modes, dissociation energies, or the KHi. The thick line shows the 
"local saturation amplitude" , the thin line the same quantity if one 
neglects the growth rate of SASI, the dotted and the dashed lines 
the contributions of the advection and the entropy stratification. 

intermediate radius between the proto-neutron star and 
the shock, approximately at r,„i„i ^ (r^h + r^)/2. This 
can be understood by the fact that higher up the shock 
parasites are efficiently stabilized by advection, while 
close to the proto-neutron star they are strongly stabi- 
lized by the entropy stratification (Fig.fTS]). The most ef- 
ficient growth of the parasites therefore takes place where 
neither advection nor stratification is strong. As the min- 
ima of the curves in Fig. [Tulare quite fiat, the parasites 
should be able to grow in a large region of the flow around 
the radius r,„ini when SASI saturates. 

The saturation amplitude predicted by our analy- 
sis of the KHi and R Ti is compared with the results 
of the simulations by [Fernandez fc ThompsonI (j2009al ) 
in Fig. [TH The RTi (thick full line) is expected to 
be the dominant parasite because it grows at smaller 
SASI amplitudes than the KHi. We note that some 
RTi structures are clear l y visib le in the simulations of 
[Fernandez fc ThompsonI (|2009aD (online edition) for e = 
0.15?;^ and e = 0-2vg, in agreement with our conclusion 
that the RTi is the dominant secondary instability. How- 
ever, these RTi structures are less obvious when e = 0. 

The amplitude in the simulations decreases by a factor 
~ 25 between e = and e = 0.2?;^, while our estimate 
decreases by a factor 15. In addition to reproducing cor- 
rectly the trend, the saturation amplitude given by this 
saturation mechanism is 15 — 50% smaller than the simu- 
lated value for all e. This is consistent with Sect. 5, where 
we found that the stability threshold of the parasites was 
~ 25 — 40% smaller than the amplitude at which their 
effect is important. Given the many approximations in- 
volved in our analytic description of the parasites, and 
the many other nonlinear effects we neglected, the com- 
parison in Fig. 1161 is considered very encouraging. 

One of our assumptions is that the background sta- 
tionary flow is unchanged, which is justifled for low sat- 
uration amplitudes but is less justified if the saturation 
amplitude is very nonlinear. The uncertainty of our an- 
alytical estimate for large saturation amplitudes is illus- 
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Fig. 16. — Saturation amplitude of SASI as a function of 
the dissociation energy t: Co mparison of the simulations by 
[Fernandez Thom pson f2009a') with the parasitic instabilities sce- 
nario. The black diamonds show the amplitu de of the saturated 
I = 1 mode in the simulations of Fer nandez fc Thompson! I l2009al ) 
(rms fluctuation of the I = 1 Legendre coefficient averaged over 
very long timcscales, after the flow has settled to a quasi-steady 
state). The empty diamonds show the saturation amplitude nor- 
malized by r^h — r*, where r^^i is the average shock radius during 
the nonlinear phase of SASI, instead of the shock radius in the sta- 
tionary flow for the black diamonds (R. Fernandez, private commu- 
nication). The thick lines are the saturation amplitude predicted 
for the most unstable mode of the KHi (gray line), and the RTi 
(black). The dotted (respectively dashed) line s hows the RTi sat- 
uration amplitude when a-gaai is set to in Eqs. I I19I I (respectively 
when one considers only the fundamental mode of SASI). 

trated in Fig. 16 by the empty diamonds, where the 
saturation amphtude is normalized using the averaged 
shock radius during the nonhnear phase of SASI instead 
of its value in the stationary flow. This new normaliza- 
tion brings the values from the simulations closer to the 
predicted ones (R. Fernandez, private communication). 

What is the dominant effect causing the dramatic de- 
crease of the SASI amplitude when e is increased? Us- 
ing Eq. ([T^ allows us to identify the contributions of 
the SASI growth rate Csasi , the relative amplitude of the 
SASI entropy wave V5sasi, the background entropy strat- 
ification V^o and the background advection velocity Vz- 

According to our linear stability analysis, the most un- 
stable SASI mode is the fundamental one if e is small, the 
first harmonics if e > 0.1 5wg and the second harmonics 
if e > 0.19u|. The comparison between the dashed line 
and thick full line in Fig.[Tni separated by about ~ 25%, 
shows that higher harmonics are more sensitive to the 
RTi than the fundamental mode (see also the thin lines 
in Fig. [Ti)) . This effect on V5sasi is partly due to the fact 
that entropy gradients in a SASI mode increase with fre- 
quency. It is striking that the sharp drop in the SASI 
amplitu de that is observed by [Fernandez fc ThompsonI 
()2009a[ ) coincides with the shift from the fundamental to 
the first radial overtone. Our analysis suggests however 
that this shift may be the cause of only a small fraction 
of the decrease, and should not be considered as a gen- 
eral feature of the saturation of SASI with dissociation. 
Indeed, repeating our analysis for other aspect ratios of 
the shock to star radius indicates that this shift can also 
happen much earlier (e.g. for rsh/r^, = 0.2) or not at all 
(e.g. for Tsh/r* = 0.6). 



The SASI amplitude required for marginal parasitic 
instability (urt = 0) can be compared to the saturation 
amplitude (ctrt = cTsasi)- From Fig. [16] (dotted and thick 
lines), we estimate that the lower growth rate of SASI 
for high dissociation energy is responsible for a ~ 25% 
decrease of the saturation amplitude of SASI. 

An important consequence of energy losses at the shock 
is a slower post-shock advection speed Vz- As the RTi 
develops more easily in a slow flow, the saturation am- 
plitude of SASI naturally decreases when dissociation is 
increased. This major effect contributes to a factor 4.5, 
evaluated by comparing the local saturation amplitude 
at the shock with and without dissociation (using a = 
in order to distinguish it from the effect of the change in 
the SASI growth rate). 

An additional factor 2 is due to the change in the flow 
profile, in particular to the decrease of the entropy strat- 
ification VSq which favors the growth of parasitic insta- 
bilities. 

7. DISCUSSION 

7.1. SASI or neutrino-driven convection? 

[Fernandez fc ThompsonI (j2009bn argued that the large 
amplitude / = 1 oscillations appearing in their numer- 
ical simulations including iron dissociation and a heat- 
ing function is due to neutrino-driven convection rather 
than SASI, since SASI is stabilized by iron dissociation 
according to [Fernandez fc ThompsonI (j2009a[) . Does iron 
dissociation at the shock really prevent SASI from grow- 
ing to large amplitudes in a realistic core collapse? In 
realistic simulations the compression factor at the shock 
can reach ~ 10, which corresponds to e = 0.14 in the 
present study and a SASI amplitude of 6% of the shock 
distance (rgh —'''*), quite smaller than without dissocia- 
tion (40%). We point out however that a significant frac- 
tion of this amplitude decrease may be an artifact of the 
parameterization, which changed the cooling function by 
a factor 127 in order to keep the ratio rsh/r^, constant. 
This parameterization has the great advantage of being 
insensitive to geometrical effects that may arise if the as- 
pect ratio between the shock and the cooling surface is 
changed. However cooling is then artificially low when 
dissociation is taken into account without heating. As a 
consequence, the resulting fiow profile may not be more 
realistic than the fiow profile without dissociation. As 
our analysis suggests that entropy gradients play an im- 
portant role in the saturation of SASI, we investigated 
the effect of keeping the cooling function constant when 
dissociation is varied, resulting in a change of the shock 
radius (from 2.5r* to 1.46r*, for < e < 0.2vg). By 
performing the same analysis as in Section 6, we then 
find that the saturation amplitude of SASI should de- 
crease significantly less than when the shock radius is 
kept constant : Ar/(rsh — r^) decreases by a factor 1.75 
only. Equivalently, Ar/rgh and Ar/r^ are decreased by 
a factor 3.5 and a factor 6 respectively (to be compared 
with a decrease of 15 when rgh is constant). Although 
the geometrical effects make a direct comparison difficult, 
the fact that all these numbers are significantly smaller 
than the former variation by a factor 15 suggests that 
the decrease of the cooling function, necessary to keep 
the shock radius constant, plays a key role in decreas- 
ing the saturation amplitude. More insight on this issue 
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may be gained by including the effect of neutrino heating 
in a parameterized manner, such that dissociation could 
be varied while both the cooling function and the shock 
radius are constant. This calculation is left for a future 
study. 

The numerical simulations bv lScheck et al.l (|2008[ ) sug- 
gest that SASI is able to grow to large amplitudes even 
in the presence of dissociation. These simulations are 
significantly more realistic than the set up studied here, 
since they include a realistic equation of state where dis- 
sociation is taken into account in a physical way, and 
a simplified treatment of neutr ino heating and cooling. 
They a lso differed from those bv lFernandez fc ThompsonI 
(|2009af) by their choice of a moving inner boundary mim- 
icking the proto neutron star contraction. In some of the 
models of this article (e.g. WOO), neutrino-driven convec- 
tion was artificially suppressed but still SASI oscillations 
could grow to non negligible amplitudes. It is however 
difficult to determine which difference between the two 
setups affects most importantly the saturation amplitude 
of SASI. 

Incidentally, it is worth noting that RTi mushrooms 
have been iden tified growing on th e SASI entropy gradi- 
ents in Fig. 7 of lScheck etldl (|2008f ) and were interpreted 
in their Sect. 6.1 as secondary convection. Although in 
that article convection was not recognized as an agent of 
SASI saturation, the fact that the RTi appears at a SASI 
amplitude close to the saturation amplitude is consistent 
with a parasitic mechanism of saturation. 

The interaction between SASI and neutrino-driven 
convection is complex and still poorly understood. Could 
neutrino-driven convection prevent the growth of SASI 
by breaking its mode structure ? One may argue 
that neutrino-driven convection does not feed upon the 
SASI mode energy, but rather converts free energy from 
the stationary gradients into vorticity. Could neutrino- 
driven convection feed SASI, either by creating vortic- 
ity which would enter t he advective-acoustic cycle, or b y 
creating sound waves ([Fernandez fc Thompso'i3l2009bD ? 
These difficult questions are beyond the scope of our 
study. 

7.2. Distinguishing RTi from KHi in the simulations 

The RTi is often characterized by finger-like or 
mushroom-like structures as in Fig. |6l while the KHi is 
characterized by vortices as in Fig.[2land[T0l However, in 
a complex fiow containing both entropy gradients, shear 
and advcction, RTi mushrooms may look like vortices 

(Fig.lUD. 

The following criterion may be more useful: the RTi 
should occur preferentially where the entropy perturba- 
tion is maximum, while the KHi occurs where the shear 
is maximum. In a sloshing mode these two maxima are 
very distinct: the entropy oscillation is maximum at the 
pole (where the shock speed is maximum) , while vorticity 
is maximum at the equator (where the inclination of the 
shock is maximu m). The parasitic s tructures visible in 
the simulations by Scheck et al. (2008) and in the movies 
published online by Fernand ez fc Thompso n (2009al) are 
more vigorous near the pole, in agreement with our anal- 
ysis. 

Furthermore, the RTi structure should grow preferen- 
tially on the half wavelength of a SASI mode where the 
entropy gradient is negative. By contrast, the KHi should 



grow on the whole exte nt of the SA S I wav elength. An 
inspection of Fig. 7 of IScheck et al.l (I2008D and of the 
movies bv lFernandez fc ThompsonI ()2009aD confirms this 
distinct feature of the RTi. 

7.3. Numerical resolution needed to resolve the parasites 

An interesting concern raised by this saturation mech- 
anism is that simulations should be able to resolve the 
parasitic instabilities properly in order to give reliable 
results on the nonlinear behavior of SASI. The RTi is a 
short wavelength instability, but as is shown in Sect. 4.3 
advection tends to stabilize the small scales and makes 
the RTi dominated by large scales. Entropy stratifica- 
tion on the other hand favors small scales. As in the 
set up studied here the RTi is found to develop where 
both stabilizing effects are important, it is hard to make 
any prediction for its dominant wavelength. Any conver- 
gence study should verify that the grid size allows for the 
growth of parasites. 

As an example, Sc heck et al.l ()2008f ) witnessed the 
growth of Rayleigh- Taylor mushrooms with a typical an- 
gular scale of Z ~ 20 — 30. They were able to cap- 
ture this small angular scale by using 360 angular zones 
for 180°. Most 2D simulations use a resolution with 
> 100 — 200 angular zones, and would probably resolve 
these scales (200 zones in ,Murphv fc Burrows! (I2008D. 
121 inlBurrows et al.l ()2006l), 128-192 inlMarek fc Jankal 
(|2009f) . and 60-120 zones in lOhnishi et al.l (|2006D ). How- 



ever 3D simulations may not be able to resolve su ch small 
scales. For example llwakami et al.l ()2008l 1200 9^ mostly 
use a resolution of 30 angular zones for 180 °, which may 
be too coarse to capture such a small scale behavior. We 
note that the saturatio n amplitude of the lo w-Z modes 
(Z = 1 - 3) in Fig. 16 of llwakami et al.l ()2008D is slightly 
smaller at "high resolution" (60 zones) than at "low res- 
olution" (30 zones). While this may be explained by a 
suppression of parasitic instabilities at low resolution, we 
cannot exclude that a different saturation process may 
take place in 3D, as discussed below. 

7.4. Effects of other physical ingredients 

The saturation mechanism described in this paper can 
be used to anticipate the effects of many other physical 
ingredients of the core collapse model (e.g. 3D versus 2D, 
the rotation rate, the magnetic field) although a detailed 
analysis is beyond the scope of this paper. 

• 3D versus 2D: 3D simulations allow for non- 
axisymmetric modes that are artificially forbidden 
in axisymmetric simulations, thus a greater num- 
ber of modes is available to the linear development 
of SASI. A single mode / = 1, to = often dom- 
inates in 2D, whereas 3 modes I = 1, m = 0,±1 
have the same growth r ate in 3D if the co llaps- 
ing co re does not rotate (iFqglizzo et al.ll2007D . Be- 
sides, Iwaka mi et al.l ()2008D found that the satu- 
rated mode amplitude is independent of to. We 
cannot exclude that nonlinear processes associated 
with the coupling between different mode, ignored 
in our analysis, ar e more important in 3D than 
in 2D. Contrary to llwakami et all ()2008D however, 
iBlondin fc Mezzacappal ( 20071 ) found that one spi- 
ral mode dominates the 3D dynamics. 
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Assuming the parasitic growth of instabihties is 
the dominant saturation mechanism, our analy- 
sis based on a hnear description of the parasites 
would predict the same saturation amplitudes of 
SASI in 2D or 3D. However the nonlinear behav- 
ior of the RTi is know n to diff e r in 3D and 2D 
(e.g. iGoncharovl (|200l . iCabotI (120061)') . and this 
may a ffect the saturation of SASI. Ilwakami et al.l 
(|2008D reported a smaller saturation amplitude of 
each individual SASI mode in 3D as compared to 
2D, although the numerical convergence of this re- 
sult should be further checked (Section 7.3). If 
confirmed, it would raise the following questions 
: is this difference in amplitude a consequence of 
the different non linear Rayleigh- Taylor behavior in 
3D? Or is this the signature of a different satura- 
tion mechanism based on the interaction of m ^ 
modes ? A more systematic parametric study, sim- 
ilar to iFernandez fc ThompsonI (|2009aD but in 3D, 
could help check the relevance of parasitic instabil- 
ities in 3D. 

• Rotation rate: lYamasaki fc Foglizzd (120081 ) have 
shown that rotation increases the growth rate of 
the spiral modes rotating in the same direction 
as the steady flow, while stabilizing the counter- 
rotating ones. If the rotation is strong enough, 
a sing le spiral mode dominates the evolution of 
SASI (jBlondin fc Mezz acappa" 2007t Ilwakami et al.1 
I2008f ). According to our analysis (Eq. ^9^), the 
larger growth rate of the spiral mode could lead to 
a larger saturation amplitude of SASI. Neverthe- 
less, a detailed calculation using the exact entropy 
and vorticity profiles of the SASI eigenmodes in a 
rotating flow is required in order to make an accu- 
rate prediction. 

• Magnetic field strength: The effect of the mag- 
netic fleld on the linear phase of SAS I is yet to 
be understood (jGuilet fc Fodizzoll2010[ ). but its ef- 
fect on parasitic instabilities can already be antic- 
ipated from the point of view of the magnetic ten- 
sion which tends to prevent motions that distort 
the magnetic field lines. This effect is stabilizing 
for the perturbations with a wave vector parallel to 
the magnetic field, but does not affect those whose 
wave vector is perpendicular. One would then ex- 
pect that the magnetic field does not change the 
maximum RTi growth rate, but selects RTi modes 
with a wave vector perpendicular to the field lines. 
In contrast, the KHi can be suppressed if the mag- 
netic field along the direction of the transverse ve- 
locity is strong enough. In a situation where the 
KHi were the dominant parasitic instability, a mag- 
netic field could potentially allow for a larger sat- 
uration amplitude. 

8. SUMMARY 

In this article we have developed for the first time a pre- 
dictive mechanism for the saturation of SASI. In this sce- 
nario the saturation happens when a parasitic instability 
is able to grow fast enough to compete with SASI. Two 
types of instabilities are of potential importance: the RTi 
growing on the entropy gradients created by SASI, and 



the KHi growing on the vorticity involved in the SASI 
mode. For each of these parasites, two stabilizing effects 
were found to be crucial: the entropy stratification in 
the stationary fiow and the advection of matter toward 
the neutron star. An estimate of the growth rates taking 
into account these effects has been obtained in Sect. 3 
and 4 for the KHi and RTi respectively. The satura- 
tion amplitude of a given SASI mode has been evaluated 
by comparing its growth rate with that of the parasitic 
instabilities. 

Using numerical simulations, we studied the effect of 
parasitic instabilities on the acoustic feedback in the sim- 
plified con text of the toy m odel introduced by iFoglizzol 
l|2009f ) and lSato et all (|2009D . This confirmed the idea of 
a threshold in amplitude above which the acoustic feed- 
back is reduced, that is determined by the ability of para- 
sitic instabilities to grow despite advection. A reasonable 
estimate of this threshold has been obtained by measur- 
ing the threshold of marginal stability, which is found 
to be 25 — 40% lower than the amplitude at which the 
feedback is decreased by 50%. 

The saturation mechanism by parasitic instabilities 
can reproduce the decrease of the SASI power with 
dissociation energy ob s erved in the simulations of 
IFernandez fc ThompsonI (|2009a[ ). Our amplitude esti- 
mate based on linear growth rates remains 15 — 50% 
lower than the saturation amplitude observed in their 
simulations, which is consistent with our simulation 
of the acoustic feedback in the toy-model. Further- 
more our analysis suggests that the RTi is the dom- 
inant secondary instability. This is consistent with 
the presence of RTi mushroom s in the s imulations of 
IFernandez fc ThompsonI (|2009aj ) as well as lScheck et al.l 
l!2008t ). 

The strong decrease of the SASI power, when both the 
dissociation energy is increased and the cooling is de- 
creased, can be traced back to 4 different effects that 
help the parasitic growth of the RTi: (i) the slower 
the advection velocity in the postshock fiow, the faster 
the propagation of the RTi against the fiow, (ii) the 
softer the negative entropy profile in the background 
fiow, the easier its destabilization by the SASI entropy 
wave, (iii) the slower the growth of SASI, the lower the 
threshold for competing parasites, (iv) the steeper the 
entropy wave in the SASI mode, the faster the RTi. 
The first two effects are the dominant o nes. We point 
out t hat other studies (Scheck et al. 2008; lOhnishi et al.l 
l2006f ). which included the effect of nuclear dissocia- 
tion, witnessed pow erful SASI oscillations, contrary to 
IFernandez fc Thomp son (20 (j9al). Th e se stud ies differed 
from that of 'Ferna ndez fc Th ompson (2009aD by the use 
of a realistic equation o f state, the inclusion of some heat- 
ing, and in the case of iScheck et al.l (|2Q08l ) the contrac- 
tion of the inner boundary. The dominant cause of the 
differing simulations has not been identified yet. Further 
investigations could shed more light on this question by 
applying the analysis described in this article to set ups 
that include one or several of these additional ingredi- 
ents (heating, realistic EOS, contraction of the innner 
boundary). 

Although the saturation mechanism proposed in 
this paper compares f a vorably with the results of 
IFernandez fc ThompsonI (|2009af ). it should be tested 
with results from other setups in order to confirm its va- 
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lidity, in particular in the case of 3D simulations. We pro- 
pose that future simulations should look for signs of the 
parasitic instabilities, and check that the angular resolu- 
tion is sufficient to resolve them. If confirmed, our results 
would open new perspectives for anticipating the effect 
on the SASI amplitude of other physical ingredients such 
as the equation of state, the heating rate, the rotation, 
and magnetic field of the progenitor star. They could 
also be useful as an input for analytical models studying 
the possible consequences of SASI, such as the mo del for 
gravit ational wave emission propos ed by .Murphy et al.l 
(2009( ) (see also lMarek et all (|2009( )'l. 
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